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Wc present an analysis of the decays B^ 7^*" (892)7 and B^ — > 7^*^(892)7 using a sample 
of about 383 million BB events collected with the BABAR detector at the PEP-II asymmetric en- 
ergy B factory. We measure the branching fractions 13{B° ->■ K*°'y) = (4.47 ± 0.10 ± 0.16) x 10"^ 
and B{B+ K*+-f) = (4.22 ± 0.14 ± 0.16) x 10"^ Wc constrain the direct CP asymmetry to be 
-0.033 < A{B -» K*j) < 0.028 and the isospin asymmetry to be 0.017 < Ao- < 0.116, where the 
limits arc determined by the 90% confidence interval and include both the statistical and systematic 
uncertainties. 
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In the Standard Model (SM), the decays B -> K*-^ [l| 
proceed dominantly through one-loop b s"f electro- 
magnetic penguin transitions. Some extensions of the 
SM predict new high-mass particles that can exist in the 
loop and alter the branching fractions from their SM pre- 
dictions. Previous measurements of the branching frac- 
tions 0, 0] are in agreement with and more precise 
than SM predictions 0j S S B 0, which suffer from 
large hadronic uncertainties. 

The time-integrated CP (A) and isospin (Ao_) asym- 
metries have smaller theoretical uncertainties flOll, and 
therefore provide more stringent tests of the SM. They 
are defined by: 



_4 = 

T{B K*-f) 



-T{B 



Ao- 



t(b" 



r(B" 



■K*%) + T{B- 



(1) 



(2) 



The 



HI, 



where the symbol F denotes the partial width. 
SM predictions for A are on the order of 1% 
while those for Aq- range from 2 to 10% [1, How- 
ever, new physics could alter these parameters signifi- 
cantly [l2, 13, J^], and thus precise measurements can 
constrain those models. In this letter, we report mea- 
surements of B{B° K*°i), B{B+ K*+-f), Ao_, 
and A. We use a data sample containing about 383 mil- 
lion BB events, corresponding to an integrated luminos- 
ity of 347 fb~^, recorded at a center-of-mass (CM) energy 
corresponding to the T{AS) mass. The data was taken 
with the BaBAR detector [la] at the PEP-II asymmet- 
ric e+e^ collider. We also make use of events simulated 
using Monte Carlo (MC) methods and a GEANT4 [16] 
detector simulation. These results supercede the previ- 
ous BaBAR measurements Q . 

B — > i^*7 decays are reconstructed in the following K* 
modes: K*"^K+Tr-, K*°~*KsTr°, K*+^K+Tr°, and 
K*~^ ^ KsTT^ ■ For each signal decay mode, the selection 
requirements described below have been optimized for 
the maximum statistical sensitivity of S/ \/S + B, where 
S and B are the rates for signal and background, re- 
spectively, where the assumed signal branching fraction 
is 4.0 X 10~^ 0. The dominant source of background 

,d, s, c) 



is continuum events (e"' 



qqij), with q 



that contain a high-energy photon from a 7r° or 77 decay 
or from an initial-state radiation (ISR) process. Back- 
grounds coming from BB events are mostly from higher- 
multiplicity 6 — ^ 57 decays, where one or more parti- 
cles have not been reconstructed, and from decays of one 
B K*^ mode that enter the signal selection of another 
mode by mis-reconstructing the K* meson. 

Photon candidates are identified as localized energy de- 
posits in the calorimeter (EMC) that are not associated 



with any charged track. The signal photon candidate is 
required to have a CM energy between 1.5 and 3.5 GeV, 
to be well-isolated and to have a shower shape consistent 
with an individual photon 17|. In order to veto photons 
from TT*^ and rj decays, we form photon pairs composed of 
the signal photon candidate and all other photon candi- 
dates in the event. We then reject signal photon candi- 
dates consistent with coming from a tt" or 77 decay based 
on a likelihood ratio that uses the energy of the partner 
photon, and the invariant mass of the pair. 

Charged particles, except those used to form Ks can- 
didates, are selected from well-reconstructed tracks that 
have at least 12 hits in the Drift Chamber (DCH), and 
are required to be consistent with coming from the e+ 
interaction region. They are identified as or tt mesons 
by the Cherenkov angle measured in the Cherenkov pho- 
ton detector (DIRC) as well as by energy loss of the track 
(dS/dx) in the silicon vertex tracker (SVT) and DCH. 
The Ks candidates are reconstructed from two oppositely 
charged tracks that come from a common vertex. In the 
K*^ — > i^s7r° {K*^ — > KsT^'^) mode, we require the invari- 
ant mass of the pair to be 0.49 < rrij^+j^- < 0.52 GeV/c^ 
(0.48 < m^+„^ < 0.52 GeV/c^) and the reconstructed 
decay length of the Ks to be at least 9.3(10) times its 
uncertainty. 

We form tt" candidates by combining two photons 
(excluding the signal photon candidate) in the event, 
each of which has an energy greater than 30 MeV in 
the laboratory frame. We require the invariant mass 



of the pair to be 0.112 < 



< 0.15 GeV/c2 and 



0.114 < m^^ < 0.15 GeV/c^ for the K*° KsTr° 
and K*~^ if+7r° modes respectively. In order to re- 
fine the tt" three-momentum vector, we perform a mass- 
constrained fit of the two photons. 

We combine the reconstructed K and tt mesons to form 
K* candidates. We require the invariant mass of the 
pair to satisfy 0.78 < ttik+t,- < 1-1 GeV/c^, 0.82 < 
rnxs-KO < 1-0 GeV/c^, 0.79 < mx+T^o < 1.0 GeV/cP, and 
0.79 < mxsTT+ < 1-0 GeV/c^. The charged track pairs 
of the K*^ ^ K'^i:^ mode are required to originate from 
a common vertex. 

The K* and high-energy photon candidates are com- 
bined to form B candidates. We define in the CM frame 
(the asterisk denotes a CM quantity) Ai? = E^- E^^^^, 
where E'^ is the energy of the B meson candidate and 
-^boam is beam energy. The beam-energy-substituted 
mass is defined as tties = \/ E^lam ~ Ps^' where p^ is the 
momentum of the B candidate. In addition, we consider 
the helicity angle 9h of the K*, defined as the angle be- 
tween the momenta one of the daughters of the K* meson 
and the B candidate in the K* rest frame. The distri- 
bution of cos Oh is sin^ 9 for signal events. Signal events 
have AE close to zero with a gaussian resolution of ap- 
proximately 50 MeV, and an tties distribution centered 
at the mass of the B meson with a gaussian resolution of 
approximately 3 MeV/c^ . We only consider candidates in 
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FIG. 1: JTigg and AE projections of the fits. The points are data, the solid line is the fit result, the dotted line is the BB 
background, and the dash-dotted line is the continuum background. The dashed line gives the total {BB and continuum) 
contribution to the background. 



the ranges -0.3 < A£; < 0.3 GeV, rriES > 5.22 GeV/c^, 
and |cos0jj| < 0.75. To eliminate badly reconstructed 
events, we apply a loose selection criterion to the vertex 
separation (and its uncertainty) along the beam axis be- 
tween the B meson candidate and the rest of the event 
(ROE). The ROE is defined as all charged tracks and 
neutral energy deposits in the calorimeter that are not 
used to reconstruct the B candidate. 

In order to reject continuum background, we combine 
13 variables into a neural network (NN). One class of 
these variables exploits the topological differences be- 
tween isotropically distributed signal events and jet-like 
continuum events by considering correlations between the 
B meson candidate and the ROE. The other class exploits 
the fact that B meson decays tend to not conserve fla- 
vor, while continuum events tend to be flavor-conser ving . 
The discriminating variables are described in Ref. [l8|. 
Each signal mode has a separately trained neural net- 
work, whose output peaks at a value of one for signal-like 
events and zero for background-like events. A selection 
is made upon the output. 

After applying all the selection criteria, there are, on 
average, ~ 1.1 B^ /B^ candidates per event in simulated 
signal events. In events with multiple candidates, we se- 
lect the candidate with the reconstructed K* mass closest 
to the nominal K* 



mass 



21|. 



We perform an unbinned extended maximum likeli- 
hood fit to extract the signal yield, constructing a sepa- 
rate fit for each mode. Since the correlations among the 
three observables (mEs, ^-B, cos 0//)^ are small, we use 
uncorrelated probability distribution functions (PDFs) 
each representing the observables to construct the likeli- 
hood function. The likelihood function is: 



M 



£ = exp - ^ i 



N 



■ M 

E 

.1=1 



where N is the number of events, M = 3 is the num- 
ber of hypotheses (signal, continuum, and BB), and rii 
is the yield of a particular hypothesis. Vi is the prod- 
uct of one-dimensional PDFS over the three dimensions 
X, and a represents the fit parameters. All types of BB 
background are included in the BB component, which is 
suppressed by the use of cos 9h- The signal toes distri- 
bution for the Kstt^ and K*+ K+tt° modes is 
described by a Crystal Ball function [l^, which has two 
tail parameters that are fixed to values obtained from 
MC. For the K*°^K+Tr- and K*+^KsTr+ modes, the 
signal Toes distribution is parameterized as a piece-wise 
function /(cc) = exp {-{x - + aL,R{x - /x)^)), 
defined to the left (L) and right (R) of /i, which is the 
peak position of the distribution. Here, (Tl,_r and ul^b. 
are the widths and measures of the tails, respectively, to 
the left and right of the peak. We constrain ctl ~ o-r, 
which is floated, and flx a^^R to values obtained from 
MC. This same function also describes the signal AE dis- 
tribution for each mode, but with different values for the 
parameters. In addition, we allow ctl and ur to float in- 
dependently. The cos 9h distribution for the signal com- 
ponent is modeled by a 2"'^ order polynomial, with all of 
its parameters floating in the fit. For the continuum hy- 
pothesis, the Toes PDF is parameterized by an ARGUS 
function [20| . with its shape parameter floating in the 
fit. The continuum AE and cosOh shapes are modeled 
by a first- or second-order polynomial with its parame- 
ters floating in the fit. Various functional forms are used 
to describe the BB background, all parameters of which 
are taken from MC simulation and held fixed. All of the 
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TABLE L The signal reconstruction efficiency e, the fitted signal yield Ns, branching fraction, B, and CP asymmetry, A, 
for each decay mode. Errors are statistical and systematic, with the exception of e and A''s, which have only systematic and 
statistical errors, respectively. 

Mode A^g B(xlO~'^) Combined g(x 10"^) A Combined A 

K+Tv- 21.8±0.8 2400.0±55.4 4.45 ± 0.10 ± 0.17 1 ,^ -0.016 ± 0.022 ± 0.007 

r 

|4.22± 0.14 ±0.16 



KsTt" 13.0±0.9 256.0±20.6 4.66 ± 0.37 ± 0.35 
K+tt" 15.3±0.8 872.7±37.6 4.38 ± 0.19 ± 0.26 I422 + OI4 + OI6 +0-040 ± 0.039 ± 0.007 
Ks7v+ 20.1±0.7 759.1±33.8 4.13 ± 0.18 ± 0.16 J ' -0.006 ± 0.041 ± 0.007 



4.47 ±0.10 ±0.16 



- 0.003 ±0.017 ±0.007 



component yields are floating. 

Figure [1] and Table U show the results of the likelihood 
fit to data. The branching fractions have been obtained 
using B{r{4S) B°B°) = 0.484 ± 0.006, B{T{4S) 
B+B-) = 0.516 ± 0.006 0. Also shown are the com- 
bined branching fractions, which have been calculated 
taking into account correlated systematic errors. 

The CP asymmetry A is measured in three modes: 
K*"^K+Tr-', K*+^K+7r", and K*+^Ksn+. In each 
of these modes, the final state of the signal B meson is 
determined by its final state daughters. The fit is ac- 
complished by performing a simultaneous fit to the two 
flavor sub-samples {K* and K ) in each mode. All shape 
parameters are assumed to be flavor independent and the 
A of each component is floated in the flt. Table [T] gives 
the individual and combined A results. 

Table HIl lists the sources of systematic uncertainty for 
the branching fractions for all four modes. The "Fit 
Model" systematic incorporates uncertainties due to im- 
perfect knowledge of the normalization and shape of the 
inclusive B — > Xgj spectra, and the choice of flxed pa- 
rameters. The "Signal PDF bias" systematic uncertainty 
characterizes any bias resulting from correlations among 
the three observables, or incorrect modeling of the sig- 
nal PDFs. The remaining sources of error on the sig- 
nal efficiency are studied using control samples in the 
data. From all of these studies, we derive signal ef- 
ficiency correction factors and associated uncertainties. 
The total corrections are 0.953, 0.897, 0.919, and 0.936 
for the K*° K+n', K*^ ~* Kstt°, K*+ K+tt°, and 
K*~^ KsTr'^ modes, respectively. The systematic error 
on A comes largely from the uncertainty in the charge 
asymmetry of the hadronic interaction of the final state 
mesons with the detector material. 

We combine the branching fractions and the ratio of 
the B+ and B° lifetime t+/to = 1.071 ± 0.009 [2l| to ob- 
tain the isospin asymmetry Aq- = 0.066 ± 0.021 ± 0.022, 
which corresponds to 0.017 < Aq- < 0.116 at the 90% 
confidence interval. We also measure A{B^ — > K*^j) = 
0.018 ± 0.028 ± 0.007. The total combined CP asymme- 
try isA^ -0.003 ±0.017 ±0.007, with a 90% confidence 
interval of -0.033 <A< 0.028. 

Figure [2] shows the relativistic P-wave Breit-Wigner 
line shape fit to the sPlot [i^l of the Kn invariant mass 



distribution of data projecting out the signal component. 
For the K*° — » Kstt'^ and K*+ — > X+tt^ modes, we con- 
volve the Breit-Wigner line shape with a Gaussian with 
a width of 10 MeV (determined from MC simulation) to 
account for detector resolution. For the K'^7t~ 
and the K*^ Kstt'^ modes, the detector resolution 
is negligible. The results are consistent with the signal 
events containing only P-wave K* mesons and no other 
Kn resonances. 




mK„(GeV/c^ mK,(GeV/c^ 



FIG. 2: Fit of a single relativistic P-wave Breit-Wigner 
line shape (solid line) to the Ktt invariant mass distribu- 
tion of the sPlot of data (points). For the K*'^ Ksiv^' and 
K*+^K+n°, the Breit-Wi gner is convolved with a Gaussian 
of width 10 MeV. 

We conclude that using a sample that is almost five 
times larger than previously used, we have signficantly 
improved upon previous measurements of the B — > K*j 
decay processes [1, 0, 01 • The measured isospin- and CP- 
asymmetries and branching fractions are consistent with 
SM expectations. 
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TABLE II: Systematic errors 


(in %) of the branching 


frac- 


tions. 










Mode 










B{T{AS) ^ B^B") / B{T{4:S) 


B+B') 1.6 


1.6 


1.6 


1.6 


BB sample size 


1.1 


1.1 


1.1 


1.1 


Tracking efficiency 


1.2 


- 


0.6 


0.8 


Particle identification 


0.6 


- 


0.6 


0.2 


Photon selection 


2.2 


2.2 


2.2 


2.2 


7r reconstruction 




3.0 


3.0 




TT and r; veto 


1.0 


1.0 


1.0 


1.0 


Ks reconstruction 




0.7 




0.7 


Neural Net efficiency 


1.5 


1.0 


1.0 


1.0 


Fit Model 


0.8 


5.6 


3.1 


1.7 


Signal PDF bias 


0.9 


2.2 


1.6 


1.4 


Sum in quadrature 


3.9 


7.5 


5.7 


4.1 
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